Part III 
DYNAMICS OF LIGHT 
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[Zquations of motion of an attracting or repelling system in vibration 
about a state of rest or motion.] 


1. Letz,y,z be the rectangular coordinates of any one point A of an sisi or repelling 
system at the time t; let x +r cos x, y +r cos B, z - r cos y be the coordinates of any other point B 
of the same system, at the same moment; and let the components of the accelerating or retarding : 
action of B on A be denoted by m cos «f (r), m cos B f (r), m cos y f (r), m being a constant factor 
which may be called the mass of the point B, & the function f (r) being positive for the case of 
attraction and negative for the case of repulsion; while the distance r between the two points is 
considered as always positive: then the differential equations of the motion of the point A, if 
we attend only to the action of B upon it, are 

d*x 
di? 


^ e meos ff (r), (1.) 


—m cos a f (r), 


d?z 
dg = mos vf (r); 


and if we take account of the actions of all the points B, C, &c., on the point A, the differential 
equations of the motion of A may be then denoted as follows, 


a= S. meosaf (r), 
d 
exa S.m cos Bf (r), (2.) 
d?z 
dp S.m cos yf (r), 


* [See Cauchy, Œuvres, 11° Série, Tome ix, pp. 390—450; 1° Série, Tome rr, pp. 111-177. Also Appendix, 
Note 11, p. 638. An abstract of this manuscript appears on page 450.] 
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the sign of summation S extending to all these other points B, C, &c., each of which has its own 
set of values of m, r, «, p, y. 

If the total number of attracting or repelling points A, B, C, &c. be n & if we attend only 
to their actions upon each other, there will be then 3n simultaneous differential equations 
of the 2™¢ order including and analogous to the 3 equations (2.); & the integrals of these 
3n equations (if known) would express the 3n coordinates of the n points of the system as 
functions of the time f, involving also, besides the n masses, 6n arbitrary constants, which 
might be determined by the initial positions of the points and their initial components of 
velocities. And for every different supposition which we might make respecting these arbitrary 
initial circumstances, we should have (in general) a different result respecting the manner of 
motion of the system, though still consistent with the laws of mutual attraction or repulsion, 
as expressed by the differential equations. 


2. Imagine then that while xz, y, z and x+rcosa, y+rcos B, z4- reos y denote as before 
the coordinates of the points A & B at the time / in one possible state of motion (or rest) of 
the system, the coordinates of the same two points at the same moment in another possible 
state of motion of the same system are denoted by x+ £, y +7, 2-- C and 

x+reosa+&+AéE, y+reosB+n+Ayn, z2+reosy+f+AZ, 
respectively; and let r+ er denote the mutual distance of these two points in this new manner 
of motion, so that 


r 4- er — V (r cos « + A£)? 4- (r cos B + A)? + (r cos y + AC?. (3.) 
' Then 
Ag An AC 
reosa+Agé  reosB--A«  rcosy-- AC DAVID. UEM eee 


Hat is EL RE (uu 20100 42.20 Sees 
Thera? "ww. rter Teron te x Ld n Ds: l+e 


will denote the 3 new cosines by which the new quantity mf (r + er) is to be multiplied, in order 
to get the three new components of accelerating or retarding action of B on A; & the differential 
equations (2.) of the old motion (or rest) of the system will be changed for the new motion, 
to the following 


? 


oe + Tel m (cosa LEES, 

dy d E 

Bi am S.m m (cos p+ a), (4) 
az a A € 

dà +o - S. m (cosy aire) 


And subtracting the equations (2.) of the old motion, from the equations (4.) of the new, 
we find the following equations: 


.c AM. OE Ds (22. -fe n) 8, eE lore) 


di? l+e ’ 

d € € 

qp - S. mcos (EES) — p ) «s IE), (5) 
d? € € 

qi 7 S meosy (AED fo) g.m BIEL), 
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which express, generally and rigorously, the laws of the differences between any two possible motions 
of any one system of attracting or repelling points: or, in other words, the laws of the dynamically 
possible deviations from any one dynamically possible manner of motion of an attracting or 
repelling system. 


3. Let us now modify these general and rigorous expressions by introducing the particular 
& approximate supposition that the length and direction of the varying line which connects any 
one moving point of the system with any other differ extremely little at any arbitrary moment t in 
the second manner of motion of the system from the length and direction of the same line at the same 


eta AL 


moment in the first manner of motion. On this hypothesis the quantities — t and e will 


be extremely small, whether r be large or small, & we may neglect their dog and products; 
so that if we denote by f’ (r) the first differential coefficient of the function f (r), the equations (5.) 
will take the simplified forms, 


= S.m cos « (rf' (r) —f (r)) e-- S. mf (r): ah. 
77 S.m cos Bf std siapa ^ (6.) 
d? 


TE = S.mcosy {rf (r) -fejet 8. mp (ry 7 


in which, by (3.), Ad x A? 
" 
«—0082—' + cos B— + cosy—>. (7.) 
Such are the approximate laws of the dynamically possible deviations from one standard 
manner of motion of an attracting or repelling system, under the condition that the mutual 
distances of the several points of that system are altered extremely little by those deviations in 
length and in direction. 


4. To illustrate by an example the foregoing general results, we may remark that they 
extend to the case of the solar system, considered as a system of points of which each attracts 
every other inversely as the square of the distance. In this case f (r) 2 1-?, and the equations 
(2.) for the standard manner of motion, which may be supposed to be the actual motion of the 

et MILA GY. g POONER Ps m e Y. And if with this standard or 


system, become —; dà zSg.——— —A aa MAGN 3 Prig 

actual motion of the solar system we compare any other motion which the same system might 
possibly have, (consistently with the present values of the masses & with the present law of 
attraction,) the deviations of the latter motion from the former are connected rigorously by 
equations of the form (5.) 


- S.m cos « ( — r=) +8. m= (r ter" 
& if these deviations be imagined such as to alter very little the length and the direction of the 
distance of any one planet from any other or from the Sun, they may then be approximately 
expressed by equations of the form (6.), namely 
di 
dt? 
in which ¢ has the value (7.). 


(r 4- er)? 


, &.; 


ibi 


= —38.mcosa.; 8. iid 


&c., 
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5. Let it next be supposed that the law of mutual action 1s such that the system is capable of 
equilibrium either exactly, or at least so very nearly that the difference shall produce no sensible 
effect in the final conclusions from this supposition: and let this possible state of equilibrium be 
taken as the standard state, to which the coordinates x, y, z refer & from which the deviations or 
displacements £, 7, £ are to be measured. Then the equations (6.) will determine approximately 
the motions of a displaced point of the system ; and «, B, y, r will become independent of the time 
and will satisfy the following equations of equilibrium, 

0=S.mcosaf(r), 0=S.mcosffi(r), O=S.mcosyf(r). (8.) 

More generally, we shall suppose that the action of any one point of the system on any other is 
so small and extends sensibly to so small a distance, and that the distribution of these points in 
space is so symmetric, as to satisfy, without sensible error, the conditions 

S. m cos a! cos B? cos y? r^ f (r) d (9.) 

S.m cos at cos B? cos y* r^ f’ (r) =0, 1 
in which h, i, 2’, i” may denote any positive integers, or zero, provided that the sum of the exponents 
i, i', i" of the cosines is an odd number. That is, we shall suppose that in the state of equilibrium, 
for every point B with mass m & with coordinates x--rcosa, y+rcos B, 2+rcosy which is 
near enough to influence sensibly the point A with coordinates z, y, z, there is another near 
point C with an equal mass m and having for coordinates 

x—reosa, y-—rceosB, z—rcosy; 
or at least that the results deduced from this supposition are not liable to sensible error. It is 
worth while to observe that when the extent of sensible action is supposed small, the approxi- 
mate equations of motion of a displaced point (6.) may be obtained from the rigorous equa- 
tions (5.) by supposing only that within this small extent of sensible action the mutual distance of 
any two points changes very little in length and in direction. 


[Conditions for plane waves.| 


6. The foregoing suppositions being admitted, we shall next inquire whether under any, and 
under what conditions the approximate equations of motion of a displaced point (6.) are compatible 
with the uniform propagation of an indefinite series of plane waves, or of displacements which succeed 
each other periodically in space and time & are the same for all the points situated at any one moment 
on any one plane parallel to a given plane;in other words, we shall inquire whether & under what 
conditions the variable displacements £, 7, & can be supposed to be periodical functions of any 
one linear function, such as ua + vy + wz + st 4- 0, of the coordinates of equilibrium zx, y, z & the 
time £, consistently with the approximate equations of motion (6.). But to simplify this inquiry, 
which is not into the most general motion possible but only into the conditions of the possibility of a 
particular motion proposed, we shall consider at present only the simplest kind of periodical 

"motion; & shall therefore assume, for the laws of displacement to be examined, the following 
expressions, 
£ —£' sin (ux 4- vy 4- wz 4- st +0), 
n=7' sin (uz 4- vy -- wz 4- st - 0), (10.) 
6 — C sin (ux 4- vy 4- wz 4- st 4- 0). 


So that the question now is, whether under any, & if so then under what conditions, the constants 
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é, n`, C, u, v, w, 8, 0 can be so chosen as to satisfy the differential equations of motion (6.) by 
the expressions of displacement (10.). If we can discover such conditions, we can then determine 
the thickness (A) ofthe wave, & the slowness (u) of its propagation ; since these will be respectively 
equal to the square roots of the expressions Loi an kr à 

u? + v? + w? 8? 

7. In tbis investigation, the operation A of differencing, by which we pass from the dis- 
placement of one to the displacement of another point of the system, affects only the coordinates 
of equilibrium zx, y, z in the expressions of displacement (10.) & not the time ¢ nor the constants 
E, n`, C, U, v, W, 8, 0, which are supposed to be the same for all the points of the system, We 
are therefore to substitute for A£ in the equations (7.) and (6.) the value 


A£-  £'sin (uz vy d-wz - w Az - v Ay 4 w Az - st-- 0) 
— é` sin (uz + vy + wz -- st -- 0) 

2 

= — 9£ sin (uz + vy + we + st-- 6) (sin ent tens 


+ £' cos (uz + vy + wz + st + 0) sin (u Ax +v Ay + w Az) CH. 


= -ap (sin 20082 t oroen teres) 
+ é` cos (ux + vy + wz + st + 0) sin (ur cos «+ vr cos B + wr cos y), 
because 
E= £' sin (ux +vy +wz + st +0) 
and 
Ax=rcos«, Ay=rcosß, Az=rcosy; (12.) 


and similar substitutions are to be made for Ay and AZ. The sign of summation S refers to the 
various values of m, r, æ, B, y corresponding to the various other points B, C, &c. of the system 
which attract or repel the point A, whose coordinates when displaced are x +é, y+7, Z+ 6; so 
that if we attend to the suppositions made in the 5th paragraph and put, for abridgment, 


2 
L= S.2m (cos a? f" (r) + sin a2 rf (r)) (sin wren tune turcos) 


2 


M = S. 2m (cos B*f' (r) + sin Bir f (r)} (sin cone trams Be wr cosy) (13. 
N — 8. 2m (cos y* f" (r) - sin y* r1 f (r)} (sin eet Pur eot] * 
and 

2 

P = 8. 2mcos Bcos y (f' (r) rf (r)} (sin urere tuner 

} 2 

Q — S. 2m cos y cos (f" (r)—r-f (r)) (un trente Ee wr cosy) (14.) 
2 

R= S.2m cos «cos B (f' (r) —r-1f (r)} (un tret re tiros 4 

HMPII i 
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418 XVII. PROPAGATION OF LIGHT IN CRYSTALS [7,8 
we shall have the three following conditions to be satisfied: 


PEL - (LE Bq + QU, 
T= — (Mn+ PL+ RE), (15.) 
Dm - (NC QE Py), 


or finally, by the supposed expressions of displacement, 
(L—s*) + Ry’ + QU —0, 
(M —s8?)n' + PP + RE =0, (16.) 
(N —s*) P + QE + Pn‘ =0. 
And these are the sought conditions for the possibility of the uniform propagation of a succession 
of plane waves with displacements of the forms (10.), in a system of attracting or repelling 
points distributed as supposed above. They do not at all restrict the constant 0 but they give 
three relations between the constants u, v, w, s and those two other constants which express the 
ratios of £', 7‘, f°. 
[Values of constants for an infinitesimally strained cubic lattice. | 
8. To simplify these conditions, we shall now introduce the additional supposition of an 
equality between the quantities m, and nearly cubical arrangement of the attracting or repelling 
points in their position of equilibrium, which differs from a cubical arrangement only by three 
small and (in general) unequal dilatations or contractions in three rectangular directions: in such 
a manner that in the expressions (13.), (14.) we may put 
rcosa=(1+a,)r,cosa,, 
rcos B — (1--5,)r, cos B,, (17.) 
r C08 y — (14 c,) r, co8 y, , 
a,, b,, c, being three small constant coefficients of dilatation or contraction & r,, «,, B,, y, being 
the quantities analogous to r, «, B, y in that exactly cubical arrangement of points from which 
the arrangement of the system, in the state of equilibrium considered lately, is now supposed 
to differ little. On this hypothesis, if we neglect all terms above the second dimension with 
respect to the small dilatations or contractions a,, b,, c,, the accurate expression 
rr, V (1--a,)*cos a? + (1+6,)? cos B? 4- (1+c,)? cos y? (18.) 
deduced from the equations (17.) will give approximately 
r —r,(1--a, cos o? +b, cos B? -- c, cos y? + 3 (a, — b,)* cos a? cos B? + 3 (b, — c, )* cos B* cos y? 
+4(c,—a,)*cosy* cosa}; (19.) 
or simply, r —r, (1 4- a, cos a? +b, cos B? +c, cos y*), (20.) 
if we neglect even the squares and products of a,,6,,c,: and in this latter order of approximation, 
the equations (17.) will also give 
cos a — cos «,, (1 +a, sin a — b, cos B? — c, 00877) 
cos B — cos B, (1+, sin 8* — c, cos y? — a, cos a2), } (21.) 
cos y — C08 y, (1 4- c, sin y? — a, cos a? — b, cos £5). | 
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On the same hypothesis, we may express the rectangular coordinates r, cos x,, r, cos B,, 
r, co8 y, of the points of the cubical arrangement as follows: 
r,cosx,—i,h, r,cosB,—k,h, r,cosy,=l,h, (22.) 
h denoting the mean interval of two eme points on the length of the side of any one of the 
elementary cubical spaces, and i, , k,, l, denoting any integer numbers positive or negative or 
null, which are however to be supposed subject to the conditions 
EP D» 0, (23.) 
in order not to include in the summation of the attractive or repulsive actions of the several 
points of the system on any one particular point the action of that point itself. And when, by 
substitution of the values (20.) and (21.), we have changed the 6 expressions (13.) & (14.) for the 
quantities L, M, N, P, Q, R to other expressions of the form 
S. F (r, cos x, ,r, cos B, ,r, co& y, , r,) (24.) 
as it is possible in every case to do, we may then transform these latter expressions to others of 
the form 
E, k, L)- v, -2,-0 FG, h, k,h, lh, Vi2+k2+2h) — F (0, 0, 0, 0), (25.) 
the sign of summation X |; 7, ^ $.—«—« indicating a summation with respect to each of the 
three independent integers i, , k,, l, from — oo to oo, and the term — F (0, 0, 0, 0) being added 
in order to allow for the condition (23.). Besides, by the supposed smallness & equality of the 
mass constants m, & by the supposed uniformity of distribution of the attracting or repelling 
points & the supposed smallness of extent of their sensible influence on each other, the expres- 
sion (25.) may be transformed without sensible error into the following triple sum 
Ec, ky 1-9, -9, -o D (AVi +k? +5) -O (0) - S. O (r,), (26.) 
if we put for abridgment 


orf le NI T, V 1— x*c08 o, , r, V 1— I disino, r,) du, dx, (27.) 


It appears also to admit of proof that in virtue of the foregoing suppositions the triple sum 
(26.) may in general be very approximately expressed by the following single sum, 

8.0 (r) - X9 (i 2): (28.) 
but the investigations of Poisson* respecting similar cases of molecular action seem to show that 
we inust distinguish this last sum È from the definite integral Í 77.0 (r,) ridr, , with which we 

0 


might otherwise be tempted to confound it. Perhaps even the reduction of the triple sum (26.) 
to the single sum (28.) may in some cases sensibly err. So that it might be necessary, for some 
laws of molecular action, to employ the less simple expression 


S. F (r,cosa,, r, cos B,, 7,cosy,, r,) 
000 9E kp 2-92 D (h Vit PI) - 0 (0)— S. 0 (r,) (29.) 
and not to admit the transformation 
S. F (r, cos, , r, cos B,, r, co8 y,, r)- Xi (o `) y: (30.) 


* [Poisson, Nouveaux Mémoires de l' Académie des Sciences, Tome vi.] 
53-2 


www.rcin.org.pl 


| 420 XVII. PROPAGATION OF LIGHT IN CRYSTALS [8, 9 


but it seems probable that such cases are rare, & that the simpler expression (30.) will usually 
be very exact. 


9. The expressions r-3f'(r)—-r-3f(r) and r-f(r) 
become, by the substitution (20.), 

7,*f' (r.) —79f (r,) 4- (r,1f" (r,) 2 3r,?f" (r,) + 37,3 (r,)} (a, cosa? +b, cos B?-- c, cos y?) (31.) 
and ro f(r, - (f (r,) 27,1 f (r,)} (a, cos a? + b, cos B? +c, cos y?); (32.) 
f" (r,) denoting the second differential coefficient of the function f (r,). We have also, by (17.), 

1? cos a? =r? cos o? (1 4- 2a,), | 
1? cos B? =r? cos B? (1+ 2b,), (33.) 
1? cos y? — 1? cos y? (1+ 2c,), 
and 1? cos B cos y — 1? cos B, cos y, (1 4- b, - c,), 
1? co8 y cos «=r? cosy, cosa, (1 +c, - a), (34.) 
1? cos a. cos Bj — 7? cos a, cos B, (1 -- a, +b,); 
and, by the same substitutions, 
(sin ur cos «+ vr cos B + wr cos r)’ A (sin ur, cosa, -- vr, cos B, +wr, cos ny 
2 2 
ts (* &, COS &, -- vr, b, cos B, +wr,c, cosy, 
2 
If then we denote by L,, M,, N,, P,, Q,, R, the values of L, M, N, P, Q, Ein the cubical 


arrangement of the points, and by 
8L SL 9L 


the coefficients of the first powers of the small dilatations or contractions a, , b, , c, in the general 
expressions of these six quantities L, &c., and if we put for abridgment 


)sin (ur, cosa, - vr, cos B, 4- wr, cos y,). (35.) 
&c. 


p =" (u cosg, +v cos B, +weosy,), (36.) 


we shall have these 6 equations 
L,— S.2m(cos of’ (r,)+sina?r7!f (r,)} sin 2, 
M,=S. 2m {cos B? f" (r,) +sin B2 r7! f (r,)) sine, | (37.) 
N, 2 S. 2m (cos y? f' (r,) - sin? r,1 f (r,)} sin al 
P,=8.2mcos B,cos y, {f'(r,)—771f (r,}sin 2, 
Q, = S. 2m cosy, cos «, { }sin?, (38.) 
R,=S.2mcos«, cos B, { }sin à, 


together with 18 others of which it is sufficient to write these 4, 
om S. 6m cos o2 (f (r,)—771f (r,)}sin? 
+ S. 2m cos a! fr, f" (r,) — 3f' (r,) -3r73 f (r,)}sin 2 


+ S. mur, cosa, {cos o f" (r,) - sina? r7 f (r,)) sin 2c; (39.) 
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= S. 2m cos B? (f' (r,)—171f (r,)) sin 2 
+ S. 2m cos a? cos B? {r, f" (r,) — 3f" (r,) + 3r71 f (r,)) sin 2 
+ S. mwr, cos B, (cos a? f' (r,) -- sin a?7,1 f (r,)} sin 2;,; (40.) 
= =S. 2m cos a? cos B, cos y, (r, f" (r,) — 3f' (r,)+3r71f (r,)) sin ? (41.) 
t S. mur, cos «, cos B, cos y, {f (r,) —7,* f (r,)) sin 2c, ; 
- P, -- S. 2m cos B3 cos y, Ír, f" (r,) — 3f' (r,) + 3r, f (r,)} sin 2 


+ S. mvr, cos B? cos y, {f (r,) 2771 f (r,)) sin 2c, : (42.) 
and the 6 quantities L, M, N, P, Q, R, in the formula of the 7th paragraph, will now be repre- 
sented as follows, SL à SL 3L 
, 8a, +0, 8b, tc, 8c, , 

10. To attain a still greater simplicity, we shall suppose further that the extent of sensible 
action of one attracting or repelling point upon another is small, yet not successively small, with 
respect to the thickness of a wave, that is, with respect to the interval in space between two succes- 
sive planes on which the attracting or repelling points of the system are similarly displaced 
from their positions of equilibrium, at any one moment of time. This interval or thickness A is 
easily seen to be expressible by the formula 


L=L,+a &c. (43.) 


». 2v 
A= (44.) 
if we put for abridgment k=Vu?+ v? + w?, (45.) 


employing the symbols w, v, w to denote the same three constants of a wave as in the equations 
(10.). We shall, therefore, swppose that the product ' 


kr =20.5, or ker, =m., (46.) 


is small, for all those distances r, or r for which the function f (r), or f (r ,), or even any one of its differ- 
ential coefficients, has any sensible value; yet not so very small but that its square and perhaps 
. higher powers may be sensible. And then, if we put 
u=ka, v=kb, w=ke, (47.) 
(so that a, b, c shall denote the cosines of the angles which one of the two outer seminormals to 
the surfaces of the wave makes with the three positive semiaxes of coordinates,) we shall be 
able to develope the square of the sine of the trinomial 
t, =4(ur, cosa, -- vr, cos B, + wr, cos y,) 

=tkr, (acosa, +b cos B, +c cos y,), (48.) 
according to the ascending powers of this small product kr,, in a converging series of which 
the three first terms are 


' KE 
sin? = 1 


2 
' 


(a cos «, -- b cos B, -- c cos y,)? 


— en (acosa, 4- b cos B, -- c cos y,)* (49.) 


kêr? 
1440 


ES (a cos «,, 4- b cos B, -- c cos y, 9: 
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and in like manner we may develope sin 2, in a series of which the first terms are 


E 
6 


sin 2, = kr, (a cosa, -- bcos B, +c cos y,) - —-" (a cos «, 4-b eos B, --ceosy,. — (50.) 


11. Weshallnow neglectthe terms which are of the fourth or higher dimensions with respect 
to kr, in the 6 expressions L,, M,, N,, P,, Q,, R,; and also the terms which are of the third 
or higher dimensions with respect to the same small quantity kr, in the 6 small remaining parts 


òL b òL òL 
Me dn ' Bb, + A 


of the 6 expressions (43.). And thus we shall have 


&o., 


L-((L 4 G)a?- (R- H)U?- (Q+ 1) ck, 
M — (M -- H)U3-- (P4 I)e? (R4- G)a) 2 (51.) 
N={(N +1) c8+(Q+@)a?+(P+H) bk, 
and 
P=2Pbck?, Q=2Qcak?, R=2Rabk*, (52.) 
if we put for abridgment ; 


seh ie oe ea 
3a, 8b, "8c; 


with 8 other equations for H, J, L, M, N, P, Q, R, and if we also put 


G= G,+a (53.) 


mr, 
2 


G,=8.—f (r,) cos a?, 


mr, 


H,- 8.5 f (r,) cos BP, (54.) 


I,- 8.7: f(r)) 008 y3, 


L,- 87r f'(r,) -f (r, jooo, 


M, =N: it }cos B1, (55.) 


mr, 
N,=8. 


)eos y1, 


mr, 


P,- 8. T^ (rf (r,) - f(r,)) cos B2 cost, 


Q, 5. T } cos y? cos o, (56.) 


R,-8.—( } cos a? cos B?, 
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8G dG 8G 
ders dide 8b, 7 Pe 50, Q^ 
8H oH oH 
AN cU Sit Sc, ”? Sa, T (57.)* 
ôI 8l 81 
iow. +N,, õa, ^ 9^ 87^ 
òL òL 
Ba, =4L, 4- S. (r,)cosa$, y 7 5.6 (ri) coma oos Hj, ic = S.A (r) cosatcosy?, 
òM 
35 = 4M, + 8.4 (r,) cos f, Sc = 8.¢(r,) cos B*cos y?, SM = 8.4(r,) cos Bt cos a?, (58.) 
N 
E edi + 8.4(r,)cosy% , e o BM (Hotel Sind (v ) cos yt cos B?, 
8c, a, òb, : t 4 
and finally 
SP oQ SR 
$a, 8b, Be.) OM eos B; 067, 
ôP òM Q oN ôR ôL 
En tS ik TE: hela as We (59.) 
òP ôN 8Q ôL R 8M 
iniu T, a 29, tg bi 3a.’ 
in which 
m wT /1. @ r.) m " , 
$= Sra (ar reset ea-mryem 02). (60.) 


[Wave velocities and directions of vibrations. | 


12. Again, by the 8th paragraph, we may transform every sum of the form 


S. (r,) cos «24 cos B2* cos y2 ( (61.) 
_ into an equivalent sum of the form 
(2, , k, ,L,) Sy (r,) 2 S. y (r,) cos «2% cos B2* cos yp, (62.) 
in which (2, , 5, , l,) is a fraction depending on the integer exponents i, , k,, l, by the following 
rule 
7 1 
($,, £,, 0,) “| cos o, gin wdw, j ib, (1 — r) dee, (63.) 
0 0 


It is easily found to be a consequence of this rule that we may change the order of the 3 ex- 
ponents 2, , k,, l, in any manner, without altering the value of the fraction (?,, &,, l,);+ and 
therefore that it is sufficient to calculate the values of these 6 different combinations, 


(1,0,0), (2,0,0), (1,1,0), (3,0,0), (2,1,0), (Ll, 1,1), (64.) 
* [Equations (57.)-(60.) follow easily from results of the type 
a (r cos x) — r, cos «,, ment cos? a,, a (r cos 8) -0, ete., 
which are a direct consequence of (17.) and (19.).] 
T [(t,,5 k,, L)- r (4, +4) p (k,, +4) r (l, +4) {27T (i +k, +l, - 3)-] 
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in order to be able to apply the transformation (62.) to all those expressions of the last 
paragraph which are of the form (61.). The 6 values (64.) are easily found to be: 


au 
(1,0,0)=4; (2,0, 0)=4; aie 1 (65.) 
(3,0,0)=4; (210)24,; (LL 1)= T05: 
If then we put for abridgment 
m 
Ry=8. Er f(r), 


R,- 8. (Sf (n) - rf (r9) (66.) 


By= 8.37 f" (r,) 3f (n) rf) 


we shall have the following results: 


G,-H,-I,- HR, (67.) 
L,-M,-N,-3H,, (68.) 
P,=Q,'= RB) = By; | (69.) 
x » B =F = 2M + 3K, (70.) 
iT. -5 =12R,+15R,, | (73.) 
ERES (74.) 
i cia ap Pe Q8) 
d -3 = 5 =2R, + BR (77.) 
T=- 5p ht Ry (78.) 


And substituting these values, in the expressions of the form (53.), we find 
G — R,(1+ 2a,)+ R, (3a, +b,  c,), 
H= Ry (14-25) - R, (3b, c, +4,), (19) 
I= R,(1+2c,)+ R, (3c, +a, +b,), 
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L=3R,(1+4a,)+3R, (ba, +b, +c,), 
M —3R, (1-- 45,) - 3R, (5b, +c, ra} 

N —3R, (1-- 4c,) - 3, (5c, +a, +b,), 

P= R,(14- 25, + 2c,)+ R, (a, 4- 3b, TRII 
Q= R, (1+ 2c, - 2a,) -- R, bete (81.) 


R= R, (1+ 2a, + 2b,) -- R, (c, + 3a, + 35,); 
and therefore 


M 4 H — R$ (1 +2b,)+ R4,(3-- 155, +c, +4a,)+ 3t, (55, +c, +a,), 
N +I= RQ(1-426,) 4- E, (8 150, +a, - b,) - 3, (5c, +a, - b,), 


R+H=R,(1+2b,)+ R, (1-- 3a, - 5b, - c,) 4- Fi, (3a, + 35, va} 


L+G= Ry (1+ 2a,)+ R, (3+ ott) HM, 16/54] 
(82.) 


PI = Ry (1+ 2c,) 4- Ry (1-- 3b, + 5c, -- a) + Ry (3b, + 3c, +a,), 
Q+ G — Ry (1-- 2a,) 4- R4, (1+ 3c, 4- 5a, +6,)+ E, (3c, 4- 3a, 4- b,), 


Led RO) IQ ie ri) 0+, +) 


(83.) 


R+G= Ry (1+ 2a,)+ R, (1+3b,+c¢, + 5a,)+ E, (35, +c, + 3a,), 
P+H= Ry (1+2b,)+ R,(1+3c,+a,+5b,)+ Ry (30, +a, + 36,). 


Hence, by (51.) and (52.), 
k*L = Ry + Ry +2R,a?+2Ry (a,a?+b,b2 +06, c?) + R (3a, +b, - c,) - 3 RS (a, +b, +c,) 
+ 4R, (3a, a? -- b, b? -- c, c?) - 2h, (62,a*? — e, b? —b,c*?); (85.) 
k-* M = Ry+ R,-- 2 R,U* - 2R, (b,b? -- c, c? -- a, a2) + R, (3b, +c, - a) - 3Es (a, +b, 4- c,) 
t 4R, (3b, D? -- c, c? -- a, a?) + 2 R, (65, b? —a,c?—c,a?); (86.) 
k-3N = Ry+ R, 4 2R c +2R, (c,c?-+a,a*+b,b%) + R, (3c, +a, +b,) - SR, (a, +b, +c,) 
t 4R, (3c, c? -- a, a? -- b, b?) -- 2 R, (6c, c? —b,a? —a,b?); (87.) 
k-*P =2R, (1+ 2b, + 2c,) bc + 2R, (a, + 3b, + 3c,) be, 
kQ =2R, (1+ 2c, + 2a,) ca+ 2R, (b, + 3c, + 8a,) ca, 
kR =2R, (1+ 2a, + 2b,)ab+ 2K, (c, + 3a, + 3b,) ab. 
And hence, finally, if we can put for abridgment 
—sk + Ry+ Ry + (2Ro+ 4R) (a, a? - b, 0? - c, c?) -- (F4, R5) (a, -5,--0,) 22e, (89.) 
the equations (16.) become 
{o+ E, (a? +a, + 4a, a?) 4- R, (62, a? — c, b? — b, c?) j 


(84.) 


(88.) 


+{R, (1+ 2a, + 2b,)+ R, (3a, + 3b, +¢,)} aby’ (90.) 
+{R, (1+ 2a, + 2c,) - E, (3a, + 2c, +b,)} act? — 0, 

{o + R, (b? +, + 4b, b?) + Ry (65, 0*— a, c? —c,a*)} m" 
+{R, (1+ 2b, + 2c,) + E, (35, +3c,+a,)}bel* (91.) 
+{R, (1+ 25, + 2a,)+ R; (3b, + 3a, +c¢,)} bag` =0, 

{o + R, (c? +c, + 4c, c?) + E, (6c, c?—b, a? —a, b?)} 2 
+ {R, (1+ 2c, + 2a,)+ R, (3c, + 3a,+b,)}caé* (92.) 
+ {R, (1+ 2c, + 25,) + E, (3c, + 3b, -- a, )) cb — 0. 

HMPH — 54 
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It is to be remembered that in these 4 last equations the quantities Ry, R1, Ra and a,, b,, c, 
are constant for any one medium depending only on the distribution of the attracting or 
repelling points of the system in their positions of equilibrium and on their law of action; 
é‘,n‘, C are also constants for any one succession of plane waves and are proportional to the three 
coordinate displacements of any one vibrating point from its position of equilibrium; a, b, c 
are the cosines of the angles of direction of a normal to such a wave; and s?k-?, with which the 
auxiliary quantity o is connected by the equation (89.), is the square of the velocity wherewith 
that wave is propagated. 


13. In the particular case when 
b,0,, (93.) 


that is, when the system of points is supposed to be, in its position of equilibrium, dilated equally 
or contracted equally in the 2 rectangular directions of x and y from the standard cubical 
arrangement, the three equations (90.), (91.), (92.) become 


(o + R, (a? +a, + 4a, a?) + R; (60,a? —a,c? —c, 02) E` 

+{R, (1+ 4a) 4- BR, (6a, 4- c,)) abn 

+ CR, (1+ 2a, + 2c,)+ R, (4a, + 3c,)} act" =0; (94.) 
{o + R, (b? -- a, + 4a, b?) + R, (6a, b? — a, c? —c,a2)) «n 

+{R, (1+ 2a, + 2c,) - R, (4a, + 3c,)] bel 

+{R,(1+4a,)+ E, (6a, +¢,)} ba =0; ur OE) 
(o + E, (c? +c, - 4c, c?) - Ra (6c, c? —a,a? —a, b?) C 

+ (R, (1+ 2a, + 2c,)+ R (4a, + 3c,)} ca£ 

-F CR, (1+ 2a, + 2c,) + R, (4a, + 3c,)}cbn‘ — 0. (96.) 

And if we further suppose that these two equal dilatations or contractions in the directions of 


x & y to vanish, so that 
a,=0, b,=0, (97.) 


(which comes to taking for the standard cubical arrangement z,, y,,2, that which differs from the 
_ position of equilibrium z, y, z only by a dilatation or contraction in the direction of the axis of z, ) 
we shall have, more simply, the 3 equations following: 


(o + R,a?— Rc, b?) € + (E, + R5c,) abn + (R,+2R,c,+3R,c,) ac? — 0, 

(o+ R,b?— R,c,a?) y +(R, + E,c,) abé + (R,+2R,c,+3R,c,) bef =0, (98.) 

{o + Ric? + Ric, (1+ 4c?) + 68,0, 2) C + (R,+2R,c,+3R,c,)c (a£ +b‘) 0; 

which give, by elimination of the ratios of £', n`, £', the following cubic equation 
0— c? -- o? E, - o?c, {R — R5 -- (4h, +7R,) e?) 
+oc, (1 — c?) R; (R4 — Ro) — oc? (1 — c?) (&R2 c? + R, Ry (1+ 16c?) + 15 R3c?) 
-e(1-ep BER, +c (1 h c?)? c? (2R, + 3R,)? Ra 
={o-—c,(1—c?) Ra} (o +c, (1— c2) Riyo + Ri+ce, c? (5E, +6R,)} 

— 9 {o —c, (1 —c?) R5) c? (1 — c2) c? (R; + R5)*: (99.) 
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and the three roots of this cubic equation will give three corresponding values of the square 
s*k~ of the velocity of propagation of a plane wave by the formula 

825-2 — Ry +R, 4- c, (R,+3R,) + 2c, c? (Ry +2R,) — 20, (100.) 
to which the formula (89.) reduces itself by the supposition (97.). With respect to the three 
directions of displacement of a vibrating point corresponding to these three velocities of 
propagation of a wave, they may be deduced from the equations (98.) by determining, for each 
value of the velocity or of the auxiliary quantity c, the ratios of £', `, f°. 


14. One root of the cubic equation (99.) is evidently 
o,=C, (1 — c?) R,; (101.) 
which is easily seen to correspond, by the equations (98.), to the following relations between 
Hse’ 
ag;+bni=9, 01-0, (102.) 
that is, to a direction of displacement or vibration which is contained upon the surface of the 
wave & is perpendicular to the axis of z, that is, to the direction of dilatation or contraction of 
the system of points when in equilibrium. And the square of the velocity of a wave when 
propagated with such vibrations is, by (100.) and (101.), 


sk = Ry-- R,+¢, (Rit Ro) + 2c, (Ry +2R,+ R,)c?; (103.) 
so that it consists of a constant part, namely, 
Ry Ri +c, (R4 4- Bo), (104.) 
which is independent of the direction of the wave, and a variable part, namely, 
2c, (E  2R, +R.) c?, (105.) 


which is proportional to the square c? of the cosine of the angle between the normal to the 
surface of the wave and the direction of dilatation or contraction. 
Another root of the same cubic (99.) is, very nearly, 
o= —c,(1—c?) R, > G06.) 
. the error being of the order of the square of the dilatation or contraction, which we have already 
supposed that we may neglect. This root corresponds to the following value of the square of 
the velocity of the wave, 


sik = Ry + R,+3c, (R,+ R5) + 2c, (Rot Ri) 67; (107.) 
which consists like the former value (103.) of a constant part, in this case 
Ro+ Ri + 36, (R4 + Ro), (108.) 


and a part proportional to the square of the cosine of the angle between the direction of dilata- 
tion or contraction & the normal to the surface of the wave, namely the part 

2c, (Ry + R,) c?. (109.) 
The vibrations, for this velocity of propagation of the wave, have their directions determined 
very nearly by the following equations, deduced from (98.) and (106.), 


6 m. 5 [1,96 
‘at abe Legi Ue RO. (110.) 
Cc 
54-72 
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They are therefore very nearly tangential to the surface of the wave, & are perpendicular to the 
vibrations which correspond to the former velocity. 
Finally, the remaining root of the cubic equation (99.) is very nearly 
gg— — E, —c,c? (b. R4, + 6R,); (111.) 
it corresponds to vibrations which are perpendicular to both the two former directions of vibra- 
tion and very nearly normal to the wave, being determined very approximately by the equations 


£5 "s _ és [1 96, t 

Sms R, (Ri Rt (112.) 

and for these last vibrations the square of the velocity of propagation of the wave is very nearly 
sgk* = Ry + 32h, +c, (R,+3R,) + 2c, (RS 4- 7 E, -- 6R,) c?, (113.) 


so that it is composed of a constant and a variable part like the squares of the two former 
velocities. With respect to the small angular deviation of the third set of vibrations from the 
normal, or of the second set of vibrations from the surface of the wave, it varies very nearly as 
the sine of twice the angle between the normal and the direction of dilatation or of contraction. 
15. Results of the same kind would have followed if we had made in (89.), (90.), (91.) and 
(92.) Gye ah d, 291 : (114.) 
that is, if we had supposed the arrangement of the attracting or repelling points of the system, 
when in equilibrium, to differ from the standard cubical arrangement by two equal dilatations 
or equal contractions in the two rectangular directions of z and y. We should then have had this 
formula for the square of the velocity of the wave, 
sk — Ry 4- R4, 4- 2a, (R,+3R,) + 2a, (Ry + 2 R4) (1 — c?) — 2o; (115.) 
the auxiliary quantity c, and the ratios of £', ņ`, £', which determine the direction of vibration, 
being found from the following equations: 
0 — (o - R,a?+a, R, (1-- 4a?) - a, R, (69? — c?) é` 
+{R, + 2a, (2R, + 3R5)] aby’ - (E, + 2a, (R,+2R,)} ack’, | 
0={c+ Rb? +a, R; (1+ 46?) +a, R, (652 —2)) n` 
+{R,+ 2a, (2R,+3R,)} ab£ +{R, + 2a, (Ry +2R,)} bel’, | 
0={o+ R,c?—a, R,(1—c?)} C - {R, + 2a, (R,-- 2 R5)) c (a£ + bn’), 
which give, by elimination, 
0— (o +a, (FE, — Ryc?) (o +a, (E, c? — R,)} (o + E, +a, (1— c?) (5R, + 6R,)} 
—9{o+a, (R, — R,c*)} a? (1. —c2) c? (R,-- R,?. (117.) 
One root of this cubic is o= —a, (E, — R,c?) (118.) 
& corresponds to the following expressions for the velocity & the vibrations, 
sik? = Ry+ R,+ 2a, (Ro+ 4R,-- 3 R4) — 2a, (Ro + 2 R,-- Ra) C, 
af; by —0, 5-0; | 
another root is, nearly, o= —a, (E,c? — R) (120.) 
& corresponds to the expressions 
sik" = Ry+ R,+ 2a, (Ry+3R,+2R,) — 2a, (RS 4- R,) c?, 
BL B fi-s RIA, ^ (121) 


(116.) 


(119.) 
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and the third root is, nearly, 


o,= — R,—a, (1— c?) (5R,+6R,), (122.) 
which gives the expressions 
sik = Ry + 3R; + 2a, (Ro +8Rh,+ 985) - 2a, (Ry+7Rk,+6R,) 2, 
$5 7s $s 14- 3a Ax. | (123.) 
a b c ^ uil 


And all these results are evidently similar in kind to those obtained in the foregoing paragraph 
and lead to analogous conclusions. 


16. In general, if we put for abridgment 
of’ -at =t, 
en - bt =n", 
a£ 4 by +cC* e e. 
| b? re c? iJ ab ^ w 
P2 paa, 


ab 2c? 125. 
n= -241 Tc P ( ) 


(124.) 


we shall get the expressions 


c 
E= —af"— by" cU, 
by means of which the system of the 3 equations (90.), (91.), (92.) may easily be transformed 
into the following: 
0— £' (o +a, R, (b?+c?)—b, R, (a? +c?) +c, (Ria? — R,b?)} 
+a(c,—a,)(R,+ Ry) (b1"— 3c¢"); (126.) 
0=7" (a +b, R, (a? 4- c?) — a, Ra (b? +c?) +c, (R10? — R,a?)} 
i +b (c, —b,) (Ry + Ra) (a£ —3cD); (127.) 
0 — c£" (o + R,+(5R,+6R,) (a, a? -- b, D? -- c, c*)) 
T 3 (E, 4- R5) a£" (a, — (a, a* +b, b? -- c, c?)) 
+3(R,+ R,) ba" (b, — (a, a? -- b, b? --c,c?)). (128.)* 
Itis easy hence to perceive that either the two quantities o and ¢", or else the three quantities 
o+ R,, £' and 7‘, must be considered as small of the order of the dilatations or contractions 
a,,6,,c,; of which two cases, the first corresponds to nearly tangential & the second case to nearly 


normal vibrations. In the first case we may suppress the products of Z", c, a,, b,, c, and then 
the equations (126.), (127.), (128.) reduce themselves to the three following: 


0 — £"' (o +a, R4 (b? -- c?) — b, R (a? -- c?) -- c, (RB,a? — R,b9)) +" (c, —a,) (R4 + R,) ab; (129.) 
0=7" (o +b, R, (a? - c?) — a, R, (D -- c2) +c, (R,0*? — Rpa?) +E (c, —b,) (R,-- R,)ab; (130.) 
0— c Ry + 3 (R, + Ry) (aa, £ "- 65, n° — (a£ by") (a,a?+b,b?+ c, e*)): (131.) 


* [The terms in (90.), (91.), (92.) not involving a,, b,, c, relate to a system of revolution, the axis of revolution 
being the direction a, b, c. This is the reason for adopting the transformation (124.), which is the simplest of its kind. 
(126.), (127.), (128.) do not follow directly from (90.), (91.), (92.). A transformation of the type of (124.) is first 
employed, e.g. (90.) multiplied by c, (91.) multiplied by a, etc.] 
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and of these the third determines a plane of vibration coinciding nearly with the surface of the 
wave, while the two first determine in general two distinct directions of vibration upon that 
plane, to each of which corresponds a certain velocity of propagation determined by the same 
two equations through the help of the two values assigned by them for the auxiliary quantity c. 
In the second case the equation (128.) reduces itself to the form 


o= — R,— (BR, + 6R,) (a, a? -- b, b? -- c, c?), (132.) 
which determines this auxiliary quantity e, and therefore also, by (89.), the velocity of 
propagation of the wave when the vibrations are nearly normal to its surface, having the 
direction determined by the equations (126.), (127.), which now reduce themselves to the 
following: 

^ b: R ^ 
g" —3ac C PW. tu c,) t, 
(133.) 


n= 2be » d A da 


It is easy to prove that this line of nearly normal vibration is perpendicular to the plane (131.) 
of the two nearly TGEA vibrations, of which the equation is, when referred to the rect- 
angular coordinates £', 7‘, £', 

0— R, (a£ + bn’ cl) +3 (R, + Ry) (ab (a, — b,) (b^ — a") 

tbc (b, —c,) (en —bf‘)+ca(c,—a,) (af‘—c&)}. (134.) 

And in this plane the two nearly tangential vibrations themselves, of which the directions are 
determined by the following quadratic equation, deduced from (129.) & (130.) by elimination 
of the auxiliary quantity ø, 

0=7"2(c, —a,)ab— ¿"2 (c, —b,)ab +n" {(c, —b,) (a2 +02) — (c, —a,) (D3--c2)), (135.) 
are easily seen to be perpendicular to each other as well as to the direction (133.) of nearly 
normal vibration;* for if we consider them as represented separately by the 2 equations 


7179,51, 7379363, (136.) 
in which, by (135.), 
_ (¢,—4,) (b? +c?) - (c, — b,) (a° +c?) uates. 
git a— (e, —a,)ab > 01g2— d (137.) 


then the two normal planes to the wave, containing the two directions (135.) or (136.) & making 
with each other an angle which in the present order of approximation must be considered as 
equal to the angle between those two directions themselves, have for their two equations in 
the rectangular coordinates £', 7‘, £', 

oni — 55; =g: se 

: 3 = F 138. 

cn, — bba = ga (c5 — a0;), 488.) 

and are therefore perpendicular to each other, because, by (137.), 
c? (1 +9192) + (b — agi) (b — ag;) = 0. (139.) 

With respect to the 2 velocities of propagation of the wave corresponding to these two rect- 
angular & nearly tangential vibrations (135.) or (136.), they may be determined (as was above 


* [Three directions perpendicular to one another are given by (90.), (91.), (92.), if the cubic in o has unequal 
roots. } 
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observed) by the formula (89.) through the help of the two values of the auxiliary quantity o 
by these 2 linear equations 
o,= —a, R, (D? +c?) +b, R, (a? -- c?) — c, (R,a?— Rb") +g, (a, —c,) (R,+ ied (140.) 
o= —G, R4 (b? -- c?) - b, R, (a? -- c?) - c, (R, a? — Rib?) +9, (a, —c,) (R4 + R5) ab; 
or, jointly, by this one quadratic, 
0 — o? -- o (ER, — R,) (a, (0? -- c?) +b, (c? -- a?) +c, (a? -- b?)) 
T (R$ 4- R2) (a,b, c? - b, c, a? + c,a, b?) 
— R, R, (a? (b? + c?)? + b? (c? + a9)? + c? (a? + 52)? + 2a, b, a*b? + 2b, c, b?c? + 2c, a,c?a?), (141.) 
which may be put under the form 
[20 + (E, — Ro) (a, (b? +c?) +b, (c? +a?) +c, (a? + 02)? 
= (E, + R5)? (a? (b? + c*)? + b? (c? + a3)? +c? (a? + 52)? + 2a, b, (ab? — c?) 
+ 2b,c, (b?c? — a?) + 2c, a, (c?a? — b?)} 
= (R, + Ry? ((c, — a, * (6? +02)? (c, — b, * (a? + 02)? + 2 (e, — a) (e, — b) (a?b? — c2) 
= (E, 4- R5)? ((a, — b,)? (c? -- a3)? + (a, — c,)? (b? -- a2)? + 2 (a, —b,) (a, —c,) (b?c?—a?)}. (142.) 


17. When any two of the three quantities a, ,b,,c, are equal to each other, the last expression 
(142.) becomes an exact square, & the quadratic equation (141.) resolves itself into two linear 
factors. For example, in the case (93.) when 


b,za,; 
the roots of the quadratic (141.) may be thus separately and rationally expressed, 


0,— —4, yide vcio iid (143.) 
o= —a, (R,— Ra) + (à, —c,) Ry (1 — c?), 
and are connected by (140.) with the following separate & rational values of g; , 2, 
a b 
n ver aca m 


which are to be substituted in the equations (136.), while the equation (131.) of the plane of 
nearly tangential vibrations becomes 


0= R,C --3(R4,-4- Ra) (a, —c,) c (a£ + b"), (145.) 
& the formula (89.) for the velocity of propagation of the wave becomes 
sk" — Ry+ R,-2-2a, (RF -3.RE,-- 3R5) +c, (R,+3R,)+2(c,—a,)c?(Ry+2R,)—20; (146.) 
so that in this case the two directions of nearly tangential vibration, with the two corresponding 
velocities of propagation of the wave, are expressed in the following manner,* 
a£, + b; — 0, 01 — 0, | (147.) 
sik-* = RQ-4- R,+ 2a, (Ro - 4E, -- 3 R5) - c, (R,4- R5) - 2 (c, — a, ) (IR 4- 2 £4 4- Ra), i 


and 
S mo Riese -a) AS | 
[^ 


s B, - (148.) 
sik — Ry - R, + 2a, (Eg 4- 3.R, -- 2R,) + 3c, (y+ Ry) - 2(c, —a,) (E 4- R4) A 


* [In the equations determining the directions of the vibrations the square of (c, —a,) is aeglected.] 


www.rcin:org.pl 


432 XVII. PROPAGATION OF LIGHT IN CRYSTALS [17,18 


In the same case (93.) the direction of nearly normal vibration & the corresponding velocity 
of the wave are determined by the following formula, deduced from equations (133.) & (132.) 
combined with (124.) & (89.): 

£& "- 63 Rı+ R; 

a B—21143(a,-0) g^, 

sgk-*= Ryo + 3K, + 2a, (Ry +8R,+9R,) +c, (R,+3R,) 
These results agree with and include those of the 14th & 15 paragraphs; and give the 

following expression for the difference of the squares of the two velocities of propagation of a 
wave, corresponding to the two directions of tangential or nearly tangential vibration, 


k-s} k= 2 (a, —c,) (R, + Ra) (1 — c?). (150.) 


(149.) 


[Fresnel’s wave surface. | 
[18.] O=07+o(R,— R,) (a, (b? -- c?) +5, (c? -- a?) +c, (a? + b?)} 
+ (R2 4- R2) (a,b, c? -- b, c, a? -- c, a, D?) (a? -- 0? +c?) 
— R, R, (a? (b? + c)? + b (c? + a2)? + 2a, b, a?b? + 2b, c, b%c? + 2c, a. c*a?); 
if c, 2 0, then 
0— c? -- c (E, — R,) (a, (b? -- c?) +b, (c? --a?)) 
+ (R? + R2)a, b, c? (a? -- 5? +c?) — R R, (a? (b? + c2)? + 2a, b, a?b? + b? (c? + a3); 


also 
2e —(—s?k-? -- RE, 4- R4 4- (E, 4- 3 ,) (a, - b,)] (a? -- 0? -- c?) 4- 2 (Ry - 2 R4) (a, a? -- b, 6); 
when a=0, b=0, c=1, then c= R,a,— Rb, or c= R,b,— Fa; 
when b=0, c=0, a=1, then o= — R45, or o= RE,b,; 
when c=0, a=0, b=1, then o= — Ra, or c= f,a,. 


In the 1% case 
Stk? — (Ro + R,)=(a, + 36,)(R,+R,) or =(3a,+6,)(R,+ Ro); 
in the 2"¢ case 
=2R,a,+ R,(5a,+36,)+3R,(a,+6,) or =2Roa, + R,(5a,+b,)+ RE, (3a, 4- b,); 
in the 374 case 
—2R,b,-- R, (3a, 4- 5b,) --3R, (a, -b,)) or  —2R$5,-- E, (a, 4- 55,) 4- E, (a, 4- 3b,). 
If these 6 values of s?k-? be called (1), (2), (3), (4), (5), (6), then 
(3)—(1)22a, (Ry +2R,+ R); (5)—(2)=2b, (R,+2R, + R,); 
(6) — (4) 2 2 (b, — a, ) (Ryo - 2 R, + Rə). 
If, then, we suppose R, + 2R,+ E, —0, we shall have 
(1) - (3) - (F-- R,) (1 —a, — 3b,), 
(2) = (5) =(Ry+ R,) (1 — 3a, —b,), 
(4) = (6) =(Ry+ R,) (1 —a, — b,). 
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Is then 
a? p? e 
#k2— (2) Pe) ek)” 
that is, 
0 — stk- — s? (a? (1) + (4) +b? (2) + (4) +c? (1) + (2)} +a? (1) (4) +b? (2) (4) -- e? (1) (2)? 
a? (1) (4) 
(Riı+ Ra)? 
b? (2) (4) 
(2, + Ro)? 
(RA gsm e (0— 4a, — 4b, + 3a + 10a, b, + 33) ; 
Sum = 1 — 4a, — 45, + 3a? + 4a, b, + 30? + 2 (a, a? -- b, b?) — 2 (a2a* — 3a b, c? + b? ?); 


— a? (1— 2a, — 45, +a? -- 4a, b, + 3b?), 


=b? (1— 4a, — 2b, + 3a? + 4a, b, - b3), 


(1) 4- (4) 
2 = 2 
a RR, 2 (1—2a, — 2b )a?, 


(2) +(4)_ 


— $2 
R,+ Rk, 


2(1—2a, —b,)b?, 


(1) + (2) 
— 2 c = c A" 2. 
c R, +R, 2 (1— 2a, — 2b,)c?; 


Sum = 2 (1 — 2a, — 2b) + 2 (a, a? -- b, b?). 


Putting 
R= — R+(Ri+ R),  R,—p (RH, o-o' (R+ R), 
then : 
Mus 12 (20 41 b.) — 2p" (a, a? 4- b, b2) — 20° 
jtm plo + (2p' + 1) (a, - b,) — 2p' (a, a* +b, 5?) — 20`. 


0— 1— 4a, — 4b, + 3a? -- 4a, b, + 3b? + 2 (a, a? -- b, b?) — 2 (a? a? — 3a,, b, c? + b2b?) 
— 2 (1— 2a, — 2b, -- a, a? -- b, b?) {1+ 2p* (a, a* +b, 5?) — (2p' + 1) (a, +b,}+ 20°} 
+ {1+ 2p" (a, a? -- b, b?) — (2p` + 1) (a, +6,) -20 j? - 2) + 2i + X5;* 
where 
29-21—241-0, 


X= —4a, — 4b, -- 2 (a, a? -- b, b?) -- 4 (a, 4- b,) — 2 (a, a? -- b, b?) — 4p' (a, a? +b, b?) 
4- 2 (2p -- 1) (a, 4- b,) — 4o* + 4p' (a, a? +b, b?) — 2 (2' -- 1) (a,+6,)+40'=0, 
126-00 (1 ST 2p") (a, +6, —a,a* —b, b?) Ta, b,c* — p` (1 =a) (a, Tb, —a,a? —b,b*?y; 
0—c?-- o (R, — R) (a, - b, —a,a? — b, b?) -- a,b, c? (R+ R,)? 
— R, R, (a, +b, — a, a? — b, b?)?; 
0— o? - o (E, — Fs) (a, (b? -- c?) +b, (c* +a?)} + (E - B3) a,b,c? (a? +b? +c?) 
— R, R, (a? (D? + c*)? + 2a, b, a?b? + D? (a? + c?)); 
& this quadratic equation does in fact agree with that lately found for o. Thus Fresnel’s law 


for the square of the normal velocity of a wave in a biaxal crystal would in fact result from the 
supposition of a rectangular and nearly cubical arrangement of attracting and repelling points, 


* [Le. absolute term, terms of lst degree, terms of 2nd degree in a,, b,, o .] 


HMPII 55 
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if we suppose the law of force & the mean interval to be such that Ry+2R,+R,=0: & the 
vibrations would then be perpendicular to the plane of polarisation, at least for uniaxal and *; 
probably for biaxal crystals. 

But if, instead, we suppose Rọ and E, to be so small with respect to R, that when multiplied 
by a, or b, they may be neglected, then the quadratic in c becomes 


0— o? — c E (a, - b, — a, a? — b, D?) + R2a,b,c?; 
& the question now is, whether this quadratic can be deduced from the developement of the 
formula 
0 a? b? ce? 
ith senda ts ee! 
which gives, accordingly, 
0 — o? — o Ra, (b? -- c*) -- b, (° +a?)}+ Ria, b,c*. 

Thus we should equally deduce Fresnel's law of velocity of a wave in a biaxal crystal, by 
supposing AK, & R, very small with respect to R,; but we should now have the vibrations 
parallel to the plane of polarisation, at least for uniaxal crystals and probably for biaxal also. 

Supposing A,— —2R,—R,, we shall deduce, very nearly, Fresnel’s law of polarisation in 
à biaxal crystal if we can prove that £', ?', {` are very nearly proportional to 

a b 
o+ Ry (a, a? --b,D2) - Ria, — Rb, o+ RS (a, a3 - b,09) + Rib, — R,a,’ 
c 
o+ E, (a, a? - b, b?) — Ra (a, - b,)' 
That is, 
0 — (e + R, (a, a? -- b, b? — b, )) (c£ — a£") +a, (R, c£ + Rak’), 
0 — (o -- E, (a, a? - b, b* — a, )) (ey — ac) +b, (Rien + Rab); 
that is [(125.)], since a£" + bn‘ + cf‘ — 0 very nearly, 
0— E {o+ R, (b, b? —b,) +a, R, (b+ c?) — aby" (E. + R5)a,, 
0—7" (e + E, (a, a? —a,) +b, R, (a? - ?)) — ab£" (R, + R,)b,; 
& these do in fact agree with the equations (129.), (130.). 
But if E, and E, be very small compared with R,, these last equations give 
f aba,R, |^ c4 E,(a?—1)a,. 
T OER OTTO ^ "abb RL”? 
also 


£ (a2—1)éE + abn" E acg" + ben" . 


i a£ LUPA CM ,, G(c&—ab), b(ey —bP) _ i 
1S o— Ra, SL aD Ty Pe 18 pus 3 "NE =0? 


ag“ $ bo" 
o— Ra, o— Rb, 
is (o—R,b,)(o—R,a,+ R,a*a,) +67), R, (o — R,a,)=02 
is o?—oR, {a, (D? -- c?) +b, (a?+c*)}+ R2a,b,c? —0? Yes. 


= 0? 


is 
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Reflexion & Refraction. 


(June 27, 1838.) (See the series of sheets begun at Kilmore near Nenagh, Sept. 27th, 1836.)* 
[19.] Components of attractive action of x + Az, y + Ay, z + Az on z, y, z .. T f(r), s f(r), 


= f(r). Suppose S. f (r) 2 0 and S.7f’ (r)=0, in order to account for the equilibrium of a point 
placed at the separating surface of two ordinary media. 
To account for the change of velocity which corresponds to ordinary refraction & yet to 

leave the density almost unaltered, suppose 

S.r (r?f' (r)+4rf (r)}’ 

S.{r?f" (r) + 4rf (r)} 

S.r[r (r*f' (r) + Arf (r)) 
S.r (r*f' (r) + 4rf (r)) | 

not very great, at least for the gases; because taking the velocity in vacuo for unity, we are to 
have 


very great; 
but 


Sa: {r+ ar)*f" (rtar)+4(r+ar)f(rt+ar)}, 30-2 .(r*f' (r) - 4rf (r), 


we J 
_ S. £r3f" (r) + 6r°f! (r) +4rf (r)} 


peal S. irf (vr) + 6n (n) + arf (r)} ii 
K =A,a+A,a7+..., A= 8.(*f (r)3 4f (3 ji cus 


However, for the gases, the extreme smallness of dilatation of the ether may be sufficient to 
account for the (probably) near proportion of this to the change of velocity. 


To account for the independence of velocity on colour in vacuo, suppose 
S. (r*f' (r) + 6r°f (7)) =0; 

S.r fref’ (r) + 6r*f (0) 

S. ir*f' (r) + 6r*f (r)} 

We are then to suppose in like manner 


S.r(r5f" (r) + FOF 9 e 
S.(rf'(r)-8Srf(tr) O — 


We have '." been led by both refraction and dispersion to suppose the quotient 
| S. r (r?if' (r) + 2i - 27231 (r))' 
S. (r?if" (r) + 2i - 2r*3f (r)} 
* [These sheets seem to have been lost. If we take equations (13.), p. 417, and replace mf (r) by To and expand, 
we get, on taking average values for Az?, Az*, ete., i.e. 47°, 17*, ete., 


or at least 
very great. 


great; 


(velocity)? = 3) S . (r*f" (r)+4rf (0) — xz S . (r*f" (r) + Gr°f (r)} + ete. 


In the other medium every distance r is supposed to be increased to r+ar.] 
- 55-2 


www.rcin.org:pl 


436 XVII. PROPAGATION OF LIGHT IN CRYSTALS [19 


& it is remarkable that this quotient becomes, when the sums for all space are changed to 
sums along one radius vector & these again to single definite integrals, 


fe (rtif* (r) * 2i 4-2 pti (r)Y de 
Í P gyatra’ (r) UA 2i4 3 pino) dr , 


in which denominator — A, (r?^*?f (r)) & numerator = Ap (r2:*3f' (r) + (20 — 1) r?*?f (r)) ; so that 
Ap. r* (rf (r))! ought on this hypothesis to be greater than Aj. 7? {rf (r)). 


(June 28th.) 

Let refracting surface be plane of wy, z being negligible for incident and reflected rays, 
but positive for refracted. Also let all the vibrations be parallel to y, so that £—0, €=0. Let 
incident vibration be 


(20 — ——- 2r —— 
171—708 (271. refracted =n} CoS (Fezi), reflected amies (ri), 


the incidence being perpendicular. When z < 0, ? 77, 4-35; when 2 » 0, 4 —35,; when 2— 0, both 
expressions hold, 25 —7 -- 93. When z « 0, x, y, z are integers i, k, l, l being negative; when z > 0, 
c, y, z are integers i, k, | multiplied by 1--«, « being a very small fraction positive or negative 
& l being positive. The components of attraction on that point, which when at rest was 0, 0, 0 


& when in motion 0, 25 cos — Set 0, are, in the directions +a, 4-9, +2, 


À > 
€ mil * 
S. SECO Lg P ef) - fte S. Tero) - fo) (n- meos. 
g HC ET), (n- ants oe [er (r) -f (0) ud (n mies ^r. 
SEEEN s" ry fes. S ef t) - f) (a-nes). 
The 15! and 37? vanish; the 224 ought to be = — x) cos 7T. We ought '*.* to have 
(S) vet [E rr o open 8 vem e r o-r), 
2 
O= (ning) 8 sim =H fq HO eqq) 8 sin E E T0 
that is, 
s £0 
i-e 8 + 8). S ero - r2 P 
o- [init nde 8 +i 8 E Gf (-f o +20), 


2mt 
+ [e is the increase in r when y is changed to y 4-:j — na Cos NH | 
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The 15 is satisfied because* 
p? S (&e.) = S (&c.) 23; 
220 2<0 
the 224 becomes 
= (uo 2M tot (y) — 2e 
-hes + s) Br ero- rm 9. 
when we introduce the lately established laws of reflection and refraction [i.e. n] +3 = u (1 — 3)]. 


Now 
S F (x, y, z, r)= S F (6 -- od, k-- ak, l+al, n-- an); 
220 1-0 


2 2 , 
+a 8 EM (T CO) TOT lt nt (y. s. ue, 

1-0 n n n 1>0 n 
because S f(n)=0, S.nf'(n)=0; & the same reason which led us to establish these 2 last 
equations seems sufficient to lead us to establish also the equation S. n? f” (n) 2 0; whereby 
the condition ending the last sentence is satisfied, even without supposing, according to the 
laws of reflected and refracted vibrations n} +n} = u (91 — nz). Thus the case of perpendicular 
incidence appears to be satisfactorily explained even without introducing the law of the vis 
viva, provided that we do employ the principle of equivalent vibrations. 


[Extension of above to dispersion in a biaxal crystal. | 
(July 9%, 1838.) 


[20.] das. Ax ek oF 28. p (E2) orl, 
" | 


r= Aba + WY Ay + Abs; òx = const. x cos T (t— V), V=out+ry+vz; 


2 
Aóx — — 28x (sin =) + useless terms [i.e. terms which vanish on summation]; 


AV =cAx+rAy+vhz; o? +T? +u =w-?*; 


5o geste OY etant, ay are 
P-8. Sye fo) s (n6. Qo div, 
then 


—(L-— o?) x+ R8y + Qóz - (M — w?) dy+ Pdz+ Rz — (N — w?) 62+ Qda+ Py; 
let plane of polarisation be zy, then we have 2 different velocities w,, w, determined as follows: 


Ist... AV-oe,Az-oj!Az, 8z-0, 8220, P=0, Q=0, R=0, 
aca-e A Ce) a) 
=19 Pass (0) AI "gu Laeta To AAAA 8. 


* [This result follows from the equations in the footnote of page 435 when we use the fact that the average 
values over a hemisphere of y*z? and 2? are r*/30 and r°/6 respectively.] 
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, 2 2 T? 2 
IDs... of = 39 {20 Ay? «(e Aa "xu sto " + (00) mata) o 


T 6703 


I8 w? = w3 sensibly? 
2 2 2 
Az=Azx,(1+a,), Ay=Ay,(1+6,), Az=Az,(1+c,), r= Gedy CER LIC ae TÉ 


&, , b,, c, are small dilatations or contractions & the lower genio refer to a cubical system; 
ps. O nap tte, g LO) vray 48, (2 ES A (a, Az? +b, Aft +c, A2?) 


aes 


^ S.r f (r,) +25 (3a, +6, +¢,) S. 2 Mey; 
18 oy Aa*Ay? pe * 25, o. (o re mI NEA G 5 (e EE Az? Ay? (0, A23 +.) 


T r 2 T, f if 
_1+2a, +26, f(r, IE 3a, 4- 3b, +c, 1 fe) h 
ai 30 8. des T, 210 Ses (; ( F. à 


' (ifA=00) w-(142a,) fas. r, f(r,) - 358. (2) 


+ (32, +36, +6 as. "(E 


ae) 
o - a}=2(a,—-),)|38. eyes (ee JE 


& *.' insensible if the value of w for the cubical system be unity. 
But we must show also that 


S. ds (Aa? -- Ay?) 4- (e I (Aa? — Ay?) 
8. [Oa (L0) a 
T f. E. 


in order to account for the existence of Fresnel's law of propagation in dispersive crystals. 


is much less than " 


8A aa S S, rSf(r Pp E "(y 
r 35 A a aah ke 
LE oy Er 
=(a,—6,) i # 
| menses] 


63 (14-2a, + 2b,) S,r*f (r) + 9(1+ 6a, +6b,+c,)S, PP) +a, e ees, PE) 


duit Badia] und ferent) | 


a, 
~ 4a, 1b; meg 


todo 
* [S,7?f(r) means S. r5 f (r,), etc.] 
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which may easily be supposed to be very small. Yet this result seems to indicate that Fresnel's 
law of propagation will hold less perfectly for the more refrangible rays. 


Consider III that velocity w which corresponds to the same plane wave as in Ist case, but 
to vibrations parallel to z & *." perpendicular to axis of crystal; that is, 


II... AV 2o, Ax—w;3!Az, S4=0, 8y—0, 
ere ostro a ns 
allan (reca. sqm sen (fy arto 


etm rief] nas iG ve] 
[63 1-742) 8,79f 9 01:92, +0, 3e) S, r5 (£) 
4 (Ba, +b, +3¢,) 8, "[ (y 


t 


ot - ode 20, e) 8, [ar D ee (2 CY) - scs s [00D + (7 (4)) |] 


which is the effect of a double refraction on the square of velocity of the wave and is sensible 
on account of the greatness of the factor 


sse C) «a. 6))]- 
(July 9th, 1838.) 


I think that [the] ordinary [ray of] violet light eoo have a slightly but sensibly different 
velocity in Iceland Spar, according as its course is parallel or perpendicular to the axis & I will 
ask Prof. Lloyd to try this experimentally.* 


T7? 
— 1890020 


(July 10%, 1838.) 


aca 8 06) 6 6))]- Ie este GGG) 
NC ttis | 
| 


wi=1+ 


» 7, (34, +b, +3¢,) RES ELT. 901 
od-14 210 S, 18902292 | | 
Let 
Q2 Q”) Q4» 
hence \ 


NP 4+02=20, ifc,=a,. 
Using Rudberg’s indices as given by Powell in Phil. Trans. 1836, Part I, page 18, I find for 
Iceland Spar by calculations of this date... (wz), = 1,6806; (wz), = 1,6860. These then I con- 


* I did ask him in presence of Poggendorf and Lord Adare in Observatory (S.E.L. room) on July 10th, 
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jecture or rather conclude from the foregoing theory (if it be right) to be the two indices of 
refraction of an ordinary ray H according as it passes through Iceland Spar in one or other of 
two different directions, namely perpendicular or parallel to axis [the axis of y]. 


(July 11%, 1838.) 
But I believe that it will be more accurate to consider the observed quantity 0,00014478.. 


Q(9-- Q9 
uw [Then] (v5), = 1,69152. I am led to think it likely that the 


value of the same index for light parallel to axis would be found to be 1,6915: the difference 
being about the 200th part of either. Malus's differences between observation and calculation 
seem to have been about as great as this, sometimes. 


as being Q{ and not as 


(July 12th, 1838.) 


The comparison of those parts of the squares of the ordinary and extraordinary velocities 
which are independent of dispersion give b, — $a, & the comparison of those parts which depend 
on dispersion give b, — 2a, ; so that, by each method, the ether in Iceland Spar ought to be about 
twice as much dilated or compressed (I suspect dilated) in the direction of the axis of the crystal 
as in a direction perpendicular thereto. 

Again, supposing 1 — Q(D: 1 — Qf? = 4a, + 3b, : 6a, +b, , we see that the second changes faster 
than the first when a, recedes from 0 and b, approaches 0, if this recession and this approach be 
d gs d(1-Qj» zda, 
nearly equal. Suppose that db, = —$da,, a,» 0, da,» 0, then dO=0, ioe eae. 
';, 20; & thus we shall explain the phenomenon observed by Mitscherlich of ii extraordinary 
index changing much faster than the ordinary with heat & increasing while the crystal tended 
to become cubical.... 


[21.] How far would my theory agree with the experiments of Malus as cited by Biot, 
Traité de Physique, Vol. 1, the plane of incidence containing the axis of the crystal? 
Let AV=cAx+rAy, v=0. 
[Introducing new notation 
c—o0o,, T= wy, "m $, 


R= 7 8, f(r), Regn. E ae "£0). 


$ 


expressions for L, M, N, R are obtained: P=0, Q=0. It will be sufficient to quote L and R.] 
L=w? {(1+ 2a,) Ro+ (3a, +b, - c,) R,+3(1+4a,) Ri +3 (5a, +b, +c,) R,} 
ra (1+ 25,) Ry+ (a, -- 3b, -- c,) E, -- (1 4- 2a, +26,) R, 4- (3a, - 3b, - c,) RB} 


n NT (1+ 4a,) R$ - (5a, +b, - c,) Ri +5(1+6a,) Ri 4- 5 (7a, +b, 4-c,) Ry} 


| 2o; 2 we 


v z {(1+ 2a, + 2b,) Ry + (3a, + 3b, +c,) Ri +3(1+4a,+ 2b,) Ri +3 (ba, 4- 3b, -- c,) Ra} 


me n OV (14-40, ) Ry + (a, + 5b, -- c,) Ry + (1+ 2a,+4b,) R1 + (3a, 4- 5b, -- c,) Ry}; 
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E —2o0,0,1(14 2a, + 2b,) R, + (3a, + 3b, - c,) Ry} 


4 3 
Sev ((14-4a, 4- 2b,) Ri 4- (ba, -- 3b, -c,) Ri}— 


3 
p 4o, 0 
Ao? 


Xo? 


{(1+ 2a,+4b,) Ri 
+ (3a, + 5b, 4- c,) R5); 
0—(L—o?)8x - R0y 2 (M — w?) dy+ Ràx — (N — w?) z; 
for one set of vibrations 8z — 0, 8y — 0, w?= N; these vibrations are perpendicular to the plane 
of zy & are tangential to wave & correspond to the extraordinary ray if the axis of the crystal 
be (as before) the axis of y; for the other set of vibrations, namely those which give the ordinary 
ray, we have 


62-0, (w?— L)(o?- M)- R?, R(8y?—932?) — (M — L)8x9y. 
[w? is then found, neglecting squares and higher powers of A-?c-?.] 
c9? — (14- a, - b) (E - R,) + (3a, + 35, - c,) (R4, 4- Ro) 

—A~*w~* {(1 + 2a, + 2b,) (FE + E) + (4a, + 45, - c,) (Ry + R5)] 

+ (a, — b,) (w2—w?) {Ro + Ry — A7? 7? (2R 4- 3.R5 + R5)); 
and this expression may, with sufficient accuracy, be written 

w= Ry+ R, 4 (3a, 4- 3b, - c, ) (Ry +R) - àw? (Ri + R5) ((4a, + 45, +¢,) + (a, — b,) (92 — 92)). 
This includes the 2 expressions 
c = Ry 4- R4 4- (3a, + 3b,+c,) (B4, + R5) - àw? (Ri + R5) (5a, + 3b, - c,), 

and w= Ry + F4 4- (3a, + 35, 4- c,) (4, + R5) - A707? (Ry + Ra) (8a, + 5b, +¢,); 
which correspond (as we have seen) to the cases of an ordinary ray perpendicular and parallel 
to the axis. 


In general the square of the ordinary velocity ought, by this theory, for Iceland Spar, to 
increase as the ray recedes from the axis, the increase being proportional to the square of the 
sine of the inclination just as for the extraordinary ray, except that this increase is much less 
rapid for the ordinary ray. And the double wave surface for Iceland Spar, and for light of any 
one colour, ought to be a system of 2 oblate spheroids of revolution touching at their common poles. 


For the extraordinary light 
w= N — os {(1+ 2a,) (Ro + Ry) + (3a, +b, +3¢,) (Ry + Ra)} 

+ wi {(1+ 2b,) (Ro + R4) + (a, + 3b, + 3c,) (E, + R,)} 

— w2 Aw? {(1 + 4a,) (Ry + Ri) + (5a, +b, + 3c,) (R1 + R3} 

— w X73o 7 (Y + 45,) (Ro + Ry) + (a, + 5b, + 3e,) (Ri + Ro)} 

= Ry + R,+ (2a, + 2b, -- 3c,) (Ry + R,)—3(a,+b,+c,)rA-*w-? (Ry + R5) 

4 (a, — b,) (92 — c0?) {R, + Ra — 22707? (Ri + R;)). 
Or, referring to axis of crystal, by making «2— «2 = 202 — 1 & observing that c, —a,, we have 
for ordinary light ... 

w= Ry+ Ri} (4a, + 3b,) (B, -- Ro) — (4a, + 5b,) Aw? (Ri + R;) 
— 2 (a, — b,) àw? (Ri + R5) w2; 
for extraordinary ... 
w2= Ry+ B+ (4a, + 3b,) (R, + Rs) — (4, + 5b,) X-3«7* (R1 + Ri) 
T 2 (a, —b,) o2 (R, + R,— 217? (E, + R5)). 
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Hence w?— w= 2 (a, — b,) o2 {R, + Ra — A737? (Ri + R5)), 
c? being square of sine of inclination to axis.* 
In the notation lately used 
OP) =(Rj+ Rj)(62,-35); P= (Rt R) (4a, +5b,); Qf» (R; + Ry) (8a, +b,); 
NO + 00 
& OP sat 


The common difference of the increasing arithmetical progression Q2, Q, O is 
2(b,—a,) (Rj, + Rə). 
Making R, + R,=1 & introducing a and b with the same meanings as Biot, we have 
a?*=1+(6a,+b,)(R,+ R;) — (8a, -- b,) (Ri + R5) àw; 
b3— 14 (4a, + 3b) (R, + Ry) — (6a, + 8b,) (Rj + By) A737; 
20 
«2 — a? — (Y xs xa] c; 
Q-—2(b,—a,) (Ej- R5; Y-2(a,—b,) (E, - R5); 


Q Q 
-b= F +i X p " V lg2—p?-. 2 ; 
Q 
p) ^ 
Thus the received constructions and rules for determining the course of the ordinary and extra- 
ordinary rays through Iceland Spar are to be modified by supposing that a and b are the two 


c — a? 2+ (0 


equatorial semidiameters of my two oblate spheroids & that CERE is the common polar 
semidiameter of both. [Hamilton finds from numerical calculations that ** Malus’ experiments 
as quoted by Biot are quite inadequate to decide a question which turns on such small 


quantities.” ] 


(July 13th, 1838.) 


The recent calculations, having been carried on nearly to che conclusion without supposing 
c, — à, , apply so far to a biaxal crystal also, and show that the section of the wave surface in 
such a crystal made by any one of the 3 principal planes is composed of 2 concentric and coaxal 
ellipses, though one of them is very nearly circular, its ellipticity being proportional to the product 
of the doubly refractive & dispersive powers or to the difference of the 2 dispersions. 

To examine other sections of the biaxal wave surface & in particular to examine whether it 
has really two optic axes now supposed and indeed known to exist, let 


L-lL190—X3*4973L]0, M= , N= , P= ,0- , R= ; 
& let AV =cAx+rAy+vz; 


* [Various terms in the apnea for w and w? have been neglected. They would however have disappeared 
in the expression for w?—w?.] 
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then L9, MX», N® have all one common part, namely 
w2{(1+2a,) Rot (3a, +b, -- c,) Ry} + e {(1+ 20,) Ry + (a, + 35, - c,) Ry} 
+ o ((1 + 2c,) Ry+(a, +b, +3c,) Ry} 
— (E, - R) (1 2 (a, 92 -- 6,02 - c,92)) — R,(1—a,—b,—c,); 
while the other parts are 
w2 (3 (1+4a,) R, -- 3 (ba, +b, -- c,) R} - X? {(1+ 2a, + 2b,) R, + (3a, 4- 3b, - c) Fs] 
+ w2{(1+ 2a, + 2c,) Ry + (3a, +b, + 3c,) Ro}, 
c {3(1+4b,) R,+3(a,+5b, -- c,) Ro} + w2{(1+ 2b, + 2c,) Ry +(a, + 3b, + 3c,) Re} 
+w2{(1+ 2a, + 2b,) R, + (3a, + 3b, +c,) Ra}, 
c? (3 (1 -- 4c,) R,+3(a,+b, + 5c,) R3) +w2{(1+ 2a, + 2c,) Ry +(3a,+b,+3c,) Ro} 
+ o {(1+ 2b, + 2c,) R, + (a, - 3b, - 3c,) Ro}. 
Of these the first is 
— (a, - b, -- c,) (R4 +R.) - 2 (R,-- Ra) (a, w} +b, e? +c, o2) - (1. — a, —5, —c,) RE, 

+ 2a, (R, -- Ra) (1+ 492) - 2 ((a, +b, - c,) (B, 4- Ra) - (1 — a, 5, —c,) Ri) os? , 
the 1*t line being common to all 3, but the 224 line being peculiar to L®. As to R®, it has been 
assigned. 

Our fundamental equations of vibration are 0=(L— w?) 8z-- Rdy+Q8z = = . These give 
wt" = o, (L8z + Rdy + Q8z) — ws (Q8z + P8y +N 8), 
c" = o, (Rx + M 8y + Pà2) — wy (Q9x + Pdy+ Nz), 
wl" — o, (Ldx+ Rdy+ Q92) + w, (R8x -- M8y + P2) +w, (Qx - Pày +N 8z), 
if éE“= w, t — wre, 7° =w,dy—w,d2z, G= wp Òt + o,9y + o,82. | 


These last expressions give 
w, 5% = E“ — w, (c, E` + em" — e, C"), 


c,8y =N" — o, (o, E + wy" — wt"), 
w,82= —w, (v, E + wyn" — oC). 
ww, E" = w, (LE + Rn‘) — wz (QE + Po") + (7, (o, L+a,R+,Q) 
+ wz, (c4 Q +w P HoN) (m, E" -- o, m — e, C"), 
ww N= o, (RE + Mq") — wy (QE + Pn’) d (— 0, (9, R+ v, M +, P) 
t oy (wQ +w P t o, N)j (o, £" + vq" — o, C"), 
ww, 6° — cx, (LE + Ry") + o, (RE - M") + w (QE + P) 
—{w,(w,L+o,R+,Q)+w,(w,R+o,M "Fo, P) +w (0, Q - o, PHN) 
X (o,£ + e," — o, C). 


Put . R=Rf'a,0,, P=P'w,w,, Q-Q'w,o; 
56-2 
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then, omitting the small terms containing £" as a factor [i.e. the vibrations are assumed to be 
nearly tangential], 
0-2 (—o? - L— Qw 4 o2 (N — D+ Q' (92 — w?) + (P^ — E) w2)} 
Tto,0,9 {R —P'-N— L4 Q' (w2 —w?) + (P' - R') 92), 
0-53 (—o* +M — P'oj +05 (N — M + P* (wi 2) + (Q' — E) o3) 
t0,0,£ {RB — QN — M + P' (02—o2)4- (Q' — RB’) o2). 

But by an investigation made lately [pp. 432—434] these formulae will give Fresnel’s laws of 
velocity d: polarisation to semiaxes a, , b,, c,, provided we can have the 4 following relations: 
C+ (1— wy) (65 — 65) = (L— Qo?) (1 — wz) + o2 (N — `w) + wzw (P* — R'), 

c? + (1 — w?) (a2 — c2) = (M — Pw?) (1. — 92) + 2 (N — Pw?) + ww? (Q` — R’), 
a,—c5— R'— P' EN -Qw -(L — Q'o$) - o (P' R), 
b2—c?= R — Q` + N— P'w?—(M — P'w?) + w2 (Q` — BR’). 

Of these 4 equations, it is remarkable that any one is a consequence of the 3 others, in- 
dependently of L, M, N, P", Q`, R`, w,, w, and in virtue of the relation w?=1—w2—w?. 
They give 

È=L +M +o (P'+ Q- R)-(Lwi+ Mwz No$ +2R'w w+ 2P‘w?w? + 2Q'w2w?), 
a=M+N+w2(Q'+R'-P'*)- (Lw?+ Mw? + Nw? +2R\w2w? + ees wi na)» 
B= N+L +03 (R+ P’—Q)—(Loi+ Mw? + Nw? + 2R'w w? + 2P* w2 w? + 2Q\w? w2 
so that it only remains to examine how nearly these 3 quantities are constant. For this examina- 
tion it will be convenient to calculate 
. Lo,-Ro,TQo, Rw,+Mw,+No,, Qw,+ Po,+No,; 
or rather 
L+ Rw +w, M+P'w?+ Rw?, N+Q'w2+ Pw. 
Any one of these will give the 2 others by transformation of the letters in the cycles a, b,c and 
T, y s. 
Confining ourselves at first to the terms independent of dispersion, we have 
c= Ry+ R,+(R,+ R3) (3a,+3b,+c,); a®=Ryo+ Rk, +(R,+ Ro) (a, + 3b, + 30,); 
b? = Rot R,+(R,+ R,) (3a,+6, + 3c,). 

How nearly would this same constancy hold when we take account of eran Eo 

."neglecting R+ R; on account of there being no dispersion in a vacuum ” 

We may then determine the velocity & polarisation of light in a biaxal at ‘As Fresnel’s 
rules, (the vibrations however being parallel to the plane of polarisation & at least very nearly 
tangential to the wave) if we employ variable semi-axes a,, b,, c, to construct the auxiliary 
ellipsoid or the wave surface; these variable semi-axes depending on the colour & on the direc- 
tion, according to the formulae 


ado Rot Rue (Ry + 8-829) (a, - 3b, -- 3c,) — aS (a,wi+b,w? +c, w?), 
B- Be 2, +( ) Ga, +b, 36) - 278225 J 
=Fy+R,+( ) (3a, +3b,+c¢,)— VERON f 
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It is remarkable that the two differences 


ato 2 (R, + R- 229) (,-4,), i-es (R +R- ES) e ty, 


are independent of the direction of the ray or wave.... 


But it is very important to observe that although we have only been led to add a common 
part to each of Fresnel’s 2 values of the square of the velocity, so that the difference of those 2 
squares is still proportional to the product of the sines of the inclinations of the wave-normal to 2 fixed 
axes A & B and thus the phenomena of the rings appears to be satisfied, yet on the other hand 
these axes seem in the present theory to be independent of the colour; which is not true for light. 
And this would seem to be a fatal objection, unless it can be removed by some modification of 
my hypothesis... 


[Here follows an attempt at reconciling the values of a,, b,, c, for topaz with the above 
results. | 
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